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Abstract 

We show how the limiting fragmentation phenomenon can arise from 
the Color Glass Condensate model of high energy QCD. We consider 
the very forward rapidity region in relativistic heavy ion collisions 
and argue that in this region, nucleus-nucleus collisions are similar 
to proton-nucleus collisions (up to shadowing corrections). We then 
use the known results for proton-nucleus cross sections to show that it 
leads to the phenomenon of limiting fragmentation in the very forward 
region of heavy ion collisions as observed at RHIC. 
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The Relativistic Heavy Ion Collider (RHIC) has opened up a new frontier 
in high energy nucleus-nucleus collisions. Many exciting and new phenomena 
have been observed which have challenged theoretical models and predictions. 
Suppression of high p t hadrons in mid rapidity, increase of baryon to pion 
ratio with p t and a large, constant anisotropy at high p t are yet to be ex- 
plained satisfactorily. In the fragmentation region (very forward rapidity), 
the PHOBOS experiment has observed the so called limiting fragmen- 
tation phenomenon |2j, shown in Fig. ([I]) for two different energies, which 
clearly shows most central (0 — 6%) charged particle multiplicities are inde- 
pendent of the center of mass energy []. In this note, we show that the Color 
Glass Condensate model of high energy nuclei can lead to a semi-quantitative 
understanding of this phenomenon. 
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Figure 1: Limiting fragmentation observed at RHIC |L|. 

It has been suggested that at high energies, due to high gluon density 
effects, a hadron or nucleus is a Color Glass Condensate and can be de- 
scribed by semi classical methods ||. This approach has been applied to 
heavy ion collisions at RHIC with some success [|]]. Proton (deutron) nu- 
cleus collisions at RHIC, scheduled to begin shortly, will greatly clarify the 
role and significance of the gluon saturation at RHIC energies |J. Here we 
show that limiting fragmentation observed at RHIC can serve as yet another 

1 The error bars shown for y/s = 200 GeV data are the average of positive and negative 
uncertainties published by PHOBOS Q. 
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indication of the importance of high gluon density effects and the Color Glass 
Condensate at RHIC energies. 

Unlike the mid rapidity region, the fragmentation region (very forward 
rapidities) in a high energy heavy ion collision, is expected to be quite similar 
to high energy proton nucleus collisions, up to shadowing corrections^. This 
is because the Quark Gluon Plasma is expected to be formed only in mid 
rapidity and will not affect particle production in the forward rapidity region. 
Also, in the fragmentation region, one can treat the target nucleus as a 
dilute system of quarks and gluons while the projectile nucleus must be 
treated as a Color Glass Condensate due to its large number of gluons. This 
is formally the same as a proton nucleus system treated in || where one 
considers scattering of quarks and gluons coming from the proton on the 
dense nucleus. For definiteness, here we focus on quark nucleus scattering 
where the cross section is given by 



da qA^ q x 1 



d 2 q t dq- d?b t (2 



7T 



5(q — p ) J d 2 r t e 



tqt-r t 



X 



±r J d 2 R t Tr c (2 - U(R t + |) - U\R t - ^ - <x dipole (r t ) 



with a similar equation for gluon scattering. This is the multiple scattering 
generalization of quark gluon scattering in pQCD and is finite as q t — > due 
to higher twist effects. To relate this to nucleus nucleus scattering in the 
very forward rapidity region, we convolute this cross section with the quark 
and gluon distributions in the target nucleus 

da AA ~* qX f do qA ^ qX 

dx q f q/A (x q ) — — — — (2) 



d?q t dq d 2 b t J d?q t dq d 2 b t 

Furthermore, since we are interested in total number of produced particles 
per unit rapidity, we will integrate over the transverse momentum q t of the 
scattered quark. We emphasize the fact that this integral is finite and can 
be done exactly. It gives 

da AA r da qA ^ qX 

dx q f q/A (x q ) (3) 



dq d 2 b t J dq d 2 b t 

2 By shadowing here, we mean any modification of the nuclear parton distributions, be 
it anti-shadowing, EMC effect etc. 
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Using p = x q \J~s/2, q — k jz and k = ■^m^e' nh , taking advantage of 
5(p~ — q~) in ([I]) to do the x q integration and including scattering of gluons 
from the nucleus, we get 



da 



AA 



drj h d?b t 



x q fq/A^q) + %g Ga{, x q 



(4) 



where x q = x g = j-^ eT>h ~ ybeam an d fq/A and Ga are the quark and gluon 
distributions of the target nucleus. Experimentally, the maximum (r] h —y beam ) 
observed at RHIC is about +2 units of rapidity pj. One can use this fact 



and that x n 



1 to show that 



zmp 



exp(— 2). We therefore have 



x f. 



x g ~ e 



-2 + Tlh-Vbeam 



(5) 



and all dependences on and z drop out (this implicitly assumes that 
hadron multiplicities are dominated by low pt pions and that z does not vary 
much with p t for low pt pions) . 

In Fig. (0) we plot da/dT]hd 2 b t from eq. (|4]) shifted for normalization and 
by the beam rapidity and compare with the most central (0 — 6%) data Jl| 
from RHIC at y/s = 200. We have used GRV98 [0] parton distributions and 
the EKS98 parameterization of nuclear shadowing ||. 
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Figure 2: Limiting fragmentation from ([|) compared to data from RHIC. 
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There are a few caveats to our results. We can not predict the overall 
normalizations, only the slope and have to normalize our results to the data 
at one reference point taken to be the target beam rapidity. It is well known 
that leading order (in a s ) calculations of cross sections suffer from large scale 
dependence. Also, the scale dependence of nuclear parton distributions, and 
in particular the gluon distribution, is very poorly known due to the limited 
Q 2 coverage of fixed target experiments. The current parameterizations of 
nuclear gluon distributions are at best an educated guess. Unfortunately, our 
results are quite sensitive to the change of scale Q 2 (the scale dependence of 
distribution functions is not written out explicitly in eq. |]). We therefore fix 
this scale by requiring that eq. ([|) give a reasonable fit to the RHIC limiting 
fragmentation data at y/s ~ 20 GeV for a couple of units of rapidity. It 
turns out that Q = 2.3 GeV works well. We then use the same scale Q in 
(fj) to predict the multiplicities at higher energies of y/s = 130 GeV and 
yfs = 200 GeV. We also show the case when Q 2 = Q 2 (y) as suggested [] 
in f|. Here y = logl/x and Q 2 (y) = Q 2 s0 exp(Ay) with Q 2 s0 = 2.0 GeV 2 
at mid rapidity and A = 0.3 {|. The choice of Q 2 = 3.0 GeV 2 leads to a 
much better agreement but is disfavored by RHIC data. Alternatively 
and if one insists on keeping Q 2 Q = 2.0 GeV 2 at mid rapidity, the choice of 
A = 0.45 leads to a good agreement with the data but this value of A = 0.45 
is too large to fit the HERA data (also, choosing a x dependent scale in 
parton distributions would seem to violate the sum rules reflecting various 
conservation laws. Nevertheless, having this as the factorization scale might 
be theoretically tempting since otherwise there is really no hard scale left 
after we integrate over all hadrons transverse momenta). We will not pursue 
this further since we are not doing a detailed quantitative study here. 

There are two principle reasons why our approach should break down as 
we get closer to the mid rapidity region. First, as one goes further away from 
the target nucleus, high gluon density effects in the target nucleus become 
important. This will show as the growth in the saturation scale of the target 
nucleus (which is of order Aq C d right at the target nucleus rapidity). To 
estimate this, we use Q 2 (Arj) = Aq CD exp (Ar/). As one goes about three 
units of rapidity away from the target nucleus, its saturation scale becomes 
appreciable (~ lGeV) and one can not describe it as a dilute system of 
partons anymore [jT0|| . 

3 During completion of this note, we learned of || which however focuses on a different 
problem. 
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Another reason why this approach should break down as one gets closer 
to mid rapidity is that the classical fields of both nuclei will become strong 
and the system will be very different from a proton nucleus collision. Also, in 
the mid rapidity region one will have to include the media effects due to the 
deconfined matter presumably produced in heavy ion collisions at RHIC. The 
media effects are presently not well understood and are beyond the scope of 
this work. 

To summarize, the underlying physics of limiting fragmentation in the 
Color Glass Condensate model is that since most particles are produced 
with transverse momenta which are below the saturation scale of the projec- 
tile nucleus (in the target nucleus reference frame), their cross sections are 
transverse momentum independent (the black disk limit). Thus the rise of 
the particle multiplicities in the very forward rapidity region (near the target 
nucleus) is due to the blackness of the projectile nucleus and the growth of 
the target nucleus parton distributions with rapidity. 
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